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LARGE INVENTORY-SERVICE OPTIMIZATION IN 
CONFIGURE-TO-ORDER SYSTEMS 

DESCRIPTION 
BACKGROUND OF THE INVENTION 

5 Field of the Invention 

The present invention generally relates to computer assisted 
manufacturing processes and, more particularly, to a model for re-engineering 
a build-to-stock operation to a configure-to-order operation centered around 
"building blocks", thereby keeping inventory only at the component level. 

1 0 Background Description 

A configure-to-order (CTO) system is a hybrid of make-to-stock and 
make-to-order operations: a set of components (subassemblies) are built to 
stock whereas the end products are assembled to order. This hybrid model is 
most suitable in an environment where the time it takes to assemble the end 

1 5 product is negligible, while the production/replenishment leadtime for each 
component is much more substantial. Personal Computer (PC) manufacturing 
is a good example of such an environment By keeping inventory at the 
component level, customer orders can be filled quickly. On the other hand, 
postponing the final assembly until order arrival provides a high level of 

20 flexibility in terms of product variety, and also achieves resource pooling in 
terms of maximizing the usage of component inventory. Therefore, the CTO 
system appears to be an ideal business process model that provides both mass 
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customization and a quick response time to order fulfillment. 

Such a hybrid model is often referred to as an assemble-to-order 
(ATO) system in the research literature. In an ATO system, usually there is a 
pre-fixed set of end-product types from which customers must choose. In 

5 contrast, a CTO system takes the ATO concept one step further, in allowing 
each customer to configure his/her own product in terms of selecting a 
personalized set of components that go into the product. Aside from checking 
that the product so configured must "make sense", there is no "menu" of 
product types that limits the customer's choice. 

1 0 PC manufacturing traditionally has been a build-to-plan (or 

build-to-forecast) process, a process that is sometimes referred to as the 
"machine-type model" (MTM) operation. There is a set of end products, or 
MTMs. Demand forecasts over a future planning horizon are generated for 
each MTM, and updated periodically for each plarining cycle, typically, a 

1 5 weekly cycle. A "materials requirements planning" (MRP) type explosion 

technique is then used to determine the requirements for the components over 
the planning horizon, based on the bill-of-materials (BOM) structure of each 
end product. Because of the random variation involved in demand forecasts, 
safety stock is usually kept for each end product, as well as at each component 

20 level, in order to meet a desirable customer service level. However, holding 
finished goods inventory for any length of time is very costly in the PC 
business, where product life cycle is measured in months and price reduction 
takes place almost every other week. 

Y. Wang, "Service Levels in Production-Inventory Networks: 

25 Bottlenecks, Tradeoffs, and Optimization", Ph.D. Dissertation, Columbia 

University, 1988, applies an asymptotic result in an optimization problem to 
minimize average inventory holding cost with a constraint on the order fill- 
rate. J. M. Swaminathan and S. R. Tayur, "Stochastic Programming Models 
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for Managing Product Variety", in Quantitative Models for Supply Chain 
Management, S. Tayur, R. Ganeshan and M. Magazine (eds.), Kluwer 
Academic Publishers, Norwell, 1999, pp. 585-624, use stochastic 
programming models to study three different strategies at the assembly stage; 

5 utilizing component commonality, postponement (the "vanilla box 

approach"), and integrating assembly task design and operations. Other related 
recent works, not necessarily in the CTO setting, include Y. Aviv and A. 
Federgruen, "The Benefits of Design for Postponement", in Quantitative 
Models for Supply Chain Management, S. Tayur, R. Ganeshan and M. 

10 Magazine (eds.), Kluwer Academic Publishers, Norwell, 1999, pp. 553-584, 
A. Garg and H. L. Lee, "Managing Product Variety: An Operations 
Perspective", in Quantitative Models for Supply Chain Management, S. Tayur, 
R. Ganeshan and M. Magazine (eds.), Kluwer Academic Publishers, Norwell, 
1999, 467-490, L. Li, "The Role of Inventory in Delivery-Time Competition", 

1 5 Management Science, 38 (1 992), 1 82-1 97, and S. Mahajan and G. J. van 

Ryzin, "Retail Inventories and Consumer Choice", in Quantitative Models for 
Supply Chain Management, S. Tayur, R. Ganeshan and M. Magazine (eds.), 
Kluwer Academic Publishers, Norwell, 1999, 491-552. 

SUMMARY OF THE INVENTION 

20 The problem solved by the present invention is to migrate a 

manufacturing process from an existing operation to a Web-based CTO 
operation where customer orders will be taken from the Internet. As the CTO 
operation will eliminate the MTM inventory, the emphasis will be shifted to 
the components, or "building blocks", which will still follow the build-to-plan 

25 scheme, due to their long leadtimes, and hence still require inventory. The 

solution involves an inventory-service trade-off of the new CTO system, with 
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resulting the performance gains, in terms of reduced inventory cost and 
increased service level. There are other benefits: There will be better forecast 
accuracy through parts commonality and risk-pooling. Customer demand is 
expected to increase, as orders will no longer be confined within a restricted 

5 set of pre-configured MTMs. 

Whereas most studies in the literature focus on certain segments of the 
supply chain, modeled as simple stand-alone queues, the present invention 
aims at modeling large-scale, end-to-end enterprise supply chains, such as 
those in the PC industry. The centerpiece of the invention is a network of 

1 0 inventory queues, which integrates inventory control with the delay-capacity 
features of queues. The network configuration is determined by the bill of 
materials (BOM) structure of the end products. The solution approach is a 
decomposition-based approximation, coupled with nonlinear optimization 
using conjugate gradient search. The model can be used not only to study the 

1 5 inventory-service tradeoff, but also to compare the effects on inventory and 
service performances through changing the network configuration; i.e., the 
supply chain structure, and to identify the best configuration. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, aspects and advantages will be better 
20 understood from the following detailed description of a preferred embodiment 
of the invention with reference to the drawings, in which: 

Figure 1 is a diagram illustrating an example of a bill-of-material 
structure of seven components used to produce three end products in an 
configure-to-order system; 
25 Figures 2A and 2B, taken together, are a flow diagram of an 

implementation of the inventory optimization algorithm according to the 
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invention; 

Figures 3A and 3B, taken together, are a block diagram showing a 
hardware platform on which the invention may be implemented; 

Figure 4 is a graph of optimal days-of-supply for selected components 
5 in non-stationary demand environment for a target service level of 95%; 

Figure 5 is a graph showing a comparison between a build-to-stock 
("as-is") and an configure-to-order ("to-be") system; 

Figure 6 is a graph showing the effect of improving forecast accuracy; 

and 

1 0 Figure 7 is a graph showing the effect of product variety on inventory. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT OF THE INVENTION 

We consider a hybrid model, by which each end product is configured 
to order from a set of components, which, in turn, are built to stock. In other 

1 5 words, no finished goods inventory is kept for any end product, whereas each 

component ("building block") has its own inventory, replenished from a 
supplier following a base-stock policy. 

Referring now to the drawings, and more particularly to Figure 1, there 
is shown a diagram illustrating an example of this hybrid model. Abstractly, 

20 there are seven components comprising the set S={ 1 ,2,...,7} which are used in 
various combinations to produce three end products comprising the set 
M={ 1,2,3}. The components are maintained in inventory, but no inventory is 
maintained of end products. Each component inventory is indexed by /, i € S, 
where S denotes the set of all components. Associated with each component is 

25 a "store", where the inventory is kept. 
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ATO Environment 

First, consider the traditional ATO environment. That is, there is a 
pre-specified set of end products (i.e., machine type models), with each type 
indexed by m E M where M denotes the set of all product types. For each 
5 product type, there is a corresponding demand stream. 

Time is discrete, indexed by /, with each time unit called a period. Let 
D m {f) denote the demand of type m in period L Each unit of type m demand 
requires a subset of components, denoted S m c 5, possibly with multiple units 
of each component. From the BOM structure of the end product, we can 
1 0 identify M t as the subset of end products that all require component T as part 

of the assembly. That is, M i = {m: i £ S m }. 

There are two kinds of leadtimes; those associated with the 
components, and those associated with the end-products: 

, / 6 5: the in-bound leadtime - the time for the supplier of 

1 5 component i to replenish to store / once an order is placed. Assume 

this leadtime is known through a given distribution. For instance, a 
normal distribution with mean and variance given. 
Z, ° ut , m 6 M\ the out-bound leadtime - the time to supply a customer 

demand of type m 9 provided there is no stockout of any component 
20 i E S m . This time includes the order processing time, the 

assembly/reconfiguration time, and the transportation time to deliver 
the order. The distribution of this leadtime is also assumed known. 

We can now express the demand at store / in period t 9 for any i 6 S and 
any t = 1, 2,..., as follows: 
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Dfi) - E DJf + LT). 

meM t 



Note this is the standard MRP type of demand-leadtime offset. We derive the 
mean and variance for D t (t) as follows: 

mm = E E ep^+omC-o. a) 

meM t f 



ValDffl = E EE[2>>!)]P[C=«] 

-EfEEt^oMC 1 ^]" 



meM t \ J 



The variance calculation above assumes independence among demand classes; 
for correlated demands, simply add to the above the covariance terms. 

CTO Environment 

The above analysis is now adapted to the CTO environment. In this 
10 environment, there is no pre-specified product menu; in principle, every order 
can require a distinct set of components. Hence, we use M to denote the set of 
product/demand families that use the same set of components 5. For instance, 
M= {low-end machines, high-end machines, servers); orM= {individuals, 
small business, corporations} . . 
1 5 Suppose D m is the demand associated with product family m in a 

certain period (for simplicity, the time argument t is omitted). Suppose each 
unit (order) in D m requires a random number of units from component / , 
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denoted as X mt , which takes on non-negative integer values. The sets Mand M t 
are still defined as before, but they now take on the following form: 

S m = S-{i:X mi =0},and M t = M- {m:X mi = 0} 

We can then translate the end-product demand into demand for each 
5 component i (per period). 

D t - E 

where X m ,{k), for £ = 1, 2,..., are independent, identically distributed (i.i.d.) 
copies of X mi . Assume the mean and the variance of X m , are known. For 
instance, starting from the joint distribution of (X m ) ieS , derived from 

1 0 empirical demand data, we can derive the marginal distributions, and then the 
mean and the variance of X mi . Applying Wald's identity (S. M. Ross, 
Stochastic Processes, John Wiley&Sons, New York, 1996), we can derive: 

E[2>,] = E EEXJE[DJ, (3) 

meM t 

Var[D,.] = E (E[i)JVar[ZJ + Var[Z)JE 2 [Zj) 

meM > I \ (4) 

= E \mxjm>y^jmj-E 2 [xjE\Dj) 

meM t 

1 5 Therefore, equations (1) and (2) can be generalized as follows: 
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mm = E eexjE epj'+omC^l < 5 ) 



Var[D,.(0] = E^jEEPif+OlPtC^] 

+ E VaiCX-jE E[D M (r + {)]P[C , =«] (6) 

- E E^/EEPj^oFtc 1 ^]' 



Base-Stock Control 

As mentioned above, suppose each store / follows a base-stock policy. 
5 Let Rtf) denote the reorder point (or, base-stock level) in period t. Express the 
reorder point as follows: 

HO := m + WMt), (7) 

where k,(t) is the desired safety factor, while u,(0 and 0,(0 can be derived (via 
queuing analysis as follows: 

10 lift) = £ E[£><»], (8) 

S=t 

aj(l) = E VartD^)], (9) 
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where if := EJX^is the expected in-bound leadtime, and E[Z),(.s)] and 

Var[AO)] follow equations (1) and (2), respectively. Note that since time is 
discrete, we shall round up any real-valued tf to the next integer to be used 

in equations (8) and (9). Also note that the reason we choose to write R,(t) as 
5 in equation (7) is because we model the demand as follows: 

A(0 = nXt) + °®'Z, (10) 

where Z denotes the standard normal variate. 

To facilitate implementation, it is often desirable to translate R,(t) into 
"days of supply" (DOS), or more precisely, periods of supply. To do so, note 
1 0 that the u,(0 part of R,(t) simply translates into periods of demand. In addition, 
we can turn the ktf)OJ{t) part of R,(f) into 

^(Qq,(0 
h(0 



periods of demand. Hence, we can express R,{f) in terms of periods of demand, 
or DOS, as follows: 



15 



DOSfl) = If 



(11) 



Note the intuitively appealing form of equation (1 1), in particular the safety- 
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stock (or rather, safety time) part, which is equal to the product of the safety 
factor and the coefficient of variation (i.e., the ratio of standard deviation to 
mean) of the component demand over the (in-bound) leadtime. 

For the inventory performance measures, the expected on-hand 
5 inventory and back-order level are as follows: 

E[/,(0] = O&HiKit)), E[£,(0] = O t (t)G(kit)\ (12) 



Here the function G( # ) is defined as: 



G(;c) := E[Z-x] + = f(z-x)$(z)dz = <{>(x)-x$(x), (13) 

with Z denoting the standard normal variate, <b and $ denoting, respectively, 
1 0 the density function and the distribution function of Z, and $ := 1 - O(x). The 
function #(•) is defined as: 

H(x) := E[x - Z] + = x + G(x) = <b(x) + x$(x). (14) 
Note that the expressions in (12) are based on the following identities: 
/,(0 = [#,(') -A(0r and B,(t) = [D,(f) - Rffl; 
1 5 along with equations (7) and (10). Furthermore, we have 

P[£> (0 ^ Rffl = P[Z ^ k,(t)] = $ (*,(0), (1 5) 
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which relates to the service requirement to be elaborated in the next section. 

Next, suppose demand is stationary, i.e., for each demand class m, 
DJt) is invariant in distribution over time. Then, the mean and the variance of 
demand per period for each component i follow the expressions in (3) and (4). 
5 Consequently, (8) and (9) reduce to the following 

U/ = fa[D}, and a] = «fVar|X> ( .]. (16) 



We can then write 

R. = ^E[D^k t $ sd[Z) ( ]; (17) 



and hence, 



10 



DOS, = 



E[D f ] 



e 



(18) 



where 0, := sd[DJ/E[D,] is the coefficient of variation of the demand per 
period for component L (Hence, 8 f / ffi is the coefficient of variation of the 

demand over the leadtime if, which is consistent with the general formula in 
equation (11).) 

1 5 Sometimes it is more appropriate to adjust the demand distribution to 

account for non-negativity. Specifically, instead of D = \i + oZ, where Z is the 
standard normal variate, we should haveD = [u + oZ] + . The adjusted mean 
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follows from equation (13) as follows: 



E[D] = oE 











t \ 


z- 


-£ 




+ = oG 


-iL 













(19) 



To derive the adjusted variance, note the following: 



oo 

E{[(Z-xYf} = j(z-x) 2 $(z)dz 



= x$ (x) + $ (x)- 2x$ (x) +x 2 $ (x) 
= $(x)-xG(x), 



where the last equation makes use of equation (13). Hence, 



Var[Z>] = 
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(20) 



For moderately large x (say, x £ 2), from equation (13), we have G(-x)= x, 
and hence 
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E[D> E[D], Var[D]« Var[D], 



from equations (19) and (20). Therefore, the above adjustment is only needed 
when the coefficient of variation of the demand, a/ji, is relatively large, say, 
0.5 or above. 

5 From equations (7), (17) and (1 8), it is clear that to identify the base- 

stock policy is tantamount to specifying the safety factor k x for each 
component inventory. In the following sections, we discuss how to set the 
safety factor values so as to achieve the best inventory-service performance as 
specified in the optimization problems below. 
10 For ease of exposition, we shall focus on stationary demand. For non- 

stationary demand, we can simply solve the optimization problems below 
period by period. 

Service Requirement 

To start with, consider the special case of each order of type m requires 
15 exactly one unit of component / 6 S m . Let a be the required service level, 
defined here as the off-shelf availability of all the components required to 
configure a unit of type m product, for any m. Let E i denote the event that 
component / is out of stock. Then, we require, for each end product mE M, 

PKs/J * 1-*- 
20 From the well-known inclusion-exclusion formula: 

E n = E my E m^E m^pE,) — , 

» t i<f i<j<k 
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we have, as an approximation, 



*S m 

There is another way to arrive at the above inequality. Suppose we 
express the service requirement as follows: 

Note that the left hand side in inequality (22) is, in fact, a lower bound of the 
no-stockout probability of the set of components in S m that is required to 
configure the end-product m; i.e., it is a lower bound of the desired off-shelf 
availability. This claim (of a lower bound) can be argued by using stochastic 

1 0 comparison techniques involving the notion of association. (Refer to, e.g., S. 
M. Ross, Stochastic Processes, 2 nd ed., Wiley, New York (1998), for 
background materials.) Intuitively, since the component inventories are driven 
by a common demand stream {DJt)} 9 and hence positively correlated, the 
chance of missing one or several components must be less than when the 

1 5 component inventories are independent, which is what is assumed by the 
product on the left hand side of inequality (22). 
Since 

n*(ty = n [!-*(*;>] * i- 

combining the above and inequality (22), we arrive at the same inequality in 
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(21). 

In the general setting of CTO, consider demand of product family m. 
Let A c S m denote a certain configuration, which occurs in this demand stream 
with probability P(^). Then the no-stockout probability, JJ, €A )» should be 
5 weighted by ?(A). Hence, the service requirement in (22) should be changed to 

a * E m n*(*£ 

** A 

= i -E^)E*(*i) 



Since 

E m = p(^>°) •= ^ < 23 > 

in the CTO environment the service requirement in (21) can be extended to the 
10 following: 

E^.l-a, (24) 

where r m follows (23). 

Note that in the CTO environment, in particular when X m > 1, the 

stockout at component / should occur more often thanO?^ due to the batch 
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size associated with each order. In particular, (15) should be modified as 
follows (omitting the time argument t): 



X mi 

ml 



a. 



« / 



which is larger than $(£.). But this gap should be insignificant, since 



where 0 .: =o/\i i is the coefficient of variation (or, forecast error); and the batch 

size of incoming orders is usually orders of magnitude smaller when compared 
against fl t , which is the mean of demand summed over all product types 
m e M t and over the leadtime. Also note that this underestimation of the 
1 0 stockout probability is compensated by the overestimation involved in (22), 
since the latter is a lower bound of the no-stockout probability. 

We can now relate the above off-shelf availability requirement to the 
standard customer service requirements expressed in terms of leadtimes, W m . 
Suppose the required service level of type m demand is: 

15 ?[W m zwJ * a, msM, (25) 

where w m 's are given data. 

Consider type m demand. We have the following two cases: 
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(i) When there is no stockout at any store ieS m - denoting the associated 
probability as Tl 0m (t), the delay is simply Z° ut , the out-bound leadtime. 

(ii) Suppose there is a stockout at one or several stores in the subset s c S m . 
Denote the associated probability as Tls m (t). Then, the delay becomes 

I, 0 " 1 + X s , where T, is the additional delay before the missing 
components in s become available. 

Hence, we can write 



In most applications, it is reasonable to assume that 
10 L? * w m and L^>w m (27) 

both hold almost surely. For instance, this is the case when the (nominal) 

outbound leadtime is nearly deterministic and shorter than what the customers 

require (this is, after all, a precondition for-running an configure-to-order 
1 5 operation); whereas the replenish leadtime for any component is substantially 

longer, hence it is virtually impossible to meet the customer service 

requirement in the case of any component stockout. 

Clearly, with the assumed inequalities in (27), the first probability on 

the right hand side of (26) becomes one, and the probabilities under the 
20 summation all become zero. Consequently, the leadtime service requirement in 

(25) reduces to the off-shelf availability % 0m , and a lower bound of the latter is 
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the left hand side of (22), as we explained earlier. 



Inventory-Service Optimization 



10 



Our objective is to minimize the maximum stockout probability among 
all end-products subject to a given inventory budget. In view of equation (21) 
and the more general service contraint in (24), the problem can be presented as 
follows: 



where c t is the unit cost of the on-hand inventory of component, U and B > 0 is 
a given parameter, the budget limit on the expected overall inventory cost. 
Recall that Ofl(k) is the expected on-hand inventory of component / (refer to 
equation(12)) and Q l follows the specification in equation (16). 

To solve the above optimization problem, we first apply a 
transformation of variables. Let 




x. := $(*,), or *, = ®~\ x )- 



(28) 



Next, denote 



g(x) := H($-\x)). 
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We can then reformulate the problem as follows: 



min P (29) 



s.t-£ V, * p, meM (30) 



E i i (3D 

165 



5 Here, the objective value p represents the upper limit on the stockout 
probability (or, "no-fill rate") over any product types. 
Several remarks are in order: 

(i) Note that #(•) is an increasing and convex function, as evident from the 
first equality in equation (14), while is decreasing and convex (convex 

1 0 when the variable takes values on the non-negative half-line). Hence, 

g( # ) is a decreasing and convex function. 

(ii) For two- end-products, m and m' 9 ifS m *cS m and r m ,^r mi for ieS m » 

then the constraint corresponding to m becomes superfluous. We 
assume in the above formulation, all such superfluous constraints have 
1 5 already been removed through preprocessing: (Note, however, the de- 

mand for such end-products (as type m 1 ) will still contribute to the 
calculation of a, in inequality (31), if m E M r ) 

(iii) The constraint in inequality (31) is always binding, since g(«) is 
decreasing as explained in (i). 

20 The Lagrangian corresponding to the above optimization problem is: 



YOR9-2000-0041 



21 



m 

meM I ieS M 



ieS 



where A m > 0 (m E M) and T| £ 0 are the Lagrangian multipliers. Hence, the 
optimal solution can be obtained from solving the following system of non- 
linear equations: 

£ r m X m - -r\cpg>(x), ieS (32) 

meM 

£ K - 1 < 33 > 

meM 

^cp^x^B (34) 



ieS 



£ V , = p, meM and X>0 * (35) 



1 0 As there is no easy way to solve the above system of non-linear 

equations, we propose a greedy algorithm, which works as follows: We first 
discretize the problem - let A be the step size. For instance, set A = 0.01% 
(one basis points), and start with x t = 0 for all L At each step, we identify the 
variable je,„ such that increasing x„ to x„ + A will yield the smallest increase in 

15 the objective value. This corresponds to identifying whether the maximum, 

over all m£M, of the left hand side of inequality (30) will either (i) remain the 
same or (ii) increase by a positive amount. 

Consider case (i) first, i.e., there is a subset, denoted A<^S, such that for 
any / 6 A, we can increase x i to x t + A, without affecting the current (5 value in 

20 the algorithm. In this case, we pick the index i*€Aas follows: 
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z* = argmax{c I .o.[g(x I .)-g(x 1 .+A)]//' jm }. (36) 



icA 



In the complement case (ii), increasing x, (by A) for any ieS will result in a 
positive increase in the P value. In this case, we still follow (36) to pick the 
index /*, replacing A by 5. That is, we choose i* such that it will result in the 
5 smallest increase in p but in the largest reduction in inventory budget. 

The above procedure is continued until the inventory budget constraint 
is satisfied; i.e., until the left hand side of inequality (3 1) is brought down to B 
or less. 

The main computation involved in each step of the above algorithm is 
10 to evaluate the increment 

Since A is a small increment, we can approximate the above difference by the 
derivative, 

c,o,[-g'fe + A/2)]-A. 
1 5 Note that from equation (10), we have 



H'(x) = -x(j>(x) + +x4>(x) = $(*)• 



And, 
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Hence, 

*(*(*)) *(*(*)) 

To summarize, the algorithm is implemented as illustrated in the flow 
diagram of Figure 2. The process begins with an initialization phase in function 
5 block 201 in which: 

For each / E S, set x { := 0. 
For each i E S 9 set r m := P(X mi > 0) 
For each m E M, set $ m := 0. 
Set p := 0. 

1 0 At this point, the process enters the outer one of three nested processing loops. 
The first step is to initialize the set of active building blocks, A, to A = {} in 
function block 202. Then the second processing loop is entered in function 
block 203 which considers each i E S, followed by entering the third processing 
loop in function block 204 which considers each end product m that uses 

1 5 component / in its bill-of-material In function block 205, set p m := p m + r m A, 
for all m such that i E S m to perform the processing in the third nested 
processing loop. Then, in function block 206, compute 6, := max OT {P J - p. A 
determination is made in decision block 207 to determine if 6, < 0. If so, then 
add / to the set of active building blocks, A:=A + {i} 9 to complete the 

20 processing in the second nested processing loop. In function block 208, 

determine if the set of active components A is non-empty; if so, set B ;= A. 
otherwise set B :=S where B is a set of component indexes. In function block 
209, find /* := arg max, 6B {-c,a/r /nr g / (x, + A/2)}, where -#'(•) follows 
equation (37). In function block 210, set* * :=x* + A to update the no- 
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stockout probability of component /*. In function block 211, compute (3 = 
max,^ J3 m , and in function block 212, check whether inequality (3 1) is 
satisfied. If yes, stop; otherwise, in function block 213 update p„ for each m E 
M t „ set $ m := (3 m + r mi A, and go to function block 202 to complete the first 
5 processing loop. 

When there is only a single product type, the above algorithm is 
guaranteed to yield the optimal solution (modulo the discretization). For 
multiple product types, to solve the problem to optimality would require, at 
each step, the examination of all subsets of S, instead of just the individual 
1 0 components i E S. For instance, if components / and;' are involved in two 

different product types, then it is possible, at some step, to increase both x t and 
Xj to x t + A and x } + A without increasing the objective value. 

Another case of interest is when each end product has (at least) one 
unique component, which we denote as im for product m. That is, = {ij is 

1 5 a singleton set. In this case, it is quite easy to check the solution returned by the 
greedy algorithm to determine whether the solution satisfies the optimality 
conditions in equations (32) to (35). To do so, first note that equation (34) is 
always satisfied. Furthermore, equation (35) is also automatically satisfied in 
this case. For if ^.^r^. < p for some end product $€ M 9 then we can always 

20 increase the value of x p that corresponds to the unique component f 4 ; and this 
will only decrease the overall inventory budget, since g(x.) is decreasing in x. . 

Hence, we only need to check equations (32) and (33). From equation (32), 
taking into account M t = {i m }> we obtain 



X = meM, 



171 „ 
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which, along with equation (33), yields 



c t o t g'(x. ) 



Vl 



What then remains is to check whether or not the X m J s and T| derived from the 
above, along with the x/s returned from the algorithm, satisfy equation (32), for 
5 / e S\{i m , m 6 A/}, the non-unique components. 

The "dual" to the above optimization problem is to minimize the 
expected total inventory cost, subject to the stockout probability of each 
product no greater than P, a given parameter. That is, 

ieS 

10 s.t£ r^p, meM (39) 

ft*,. 



The greedy heuristic described earlier still applies. The only difference is the 
stopping criterion. Here, we stop when some (if not all) of the constraints in 
(39) become binding such that no x, values can be further increased. (That is, 
the set of end products corresponds to these binding constraints spans over the 

1 5 entire set of components S.) 

Figure 3 illustrates in block diagram form a hardware platform on 
which the algorithm illustrated in the flow diagram of Figure 2 may be 
implemented. There are illustrated two databases 301 and 302 which provide 
inputs to processor 303. Database 301 contains optimization control 

20 parameters; specifically, budget limit on expected overall inventory cost and 
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step size for the greedy algorithm. Database 302 contains manufacturing 
system parameters. These are (i) demand for end product mGM in period 
t = 1 , 2,..., T, (ii) inbound leadtime distribution for component ieS, (iii) 
outbound leadtime distribution for end product mGM, (iv) the distribution of 

5 the number of units of component i used in the assembly of end product mEM, 
(v) unit cost of component ieS 9 and (5) bill-of-materials for end product mEM. 
The processor 303 accesses these databases and first performs preprocessing 
steps of propagating demand using equations (5) and (6) and computing the 
mean and variance of the demand over the leadtime using equation (16). After 

1 0 the preprocessing, the processor performs optimization processing. This 

includes applying variable transformation using equation (28), reformulating 
the optimization problem using equations (29) to (31), and solving the 
optimization problem through the greedy algorithm as illustrated in the flow 
diagram of Figure 2. Finally, the processor enters the postprocessing phase 

15 where the reorder points are computed using equations (5), (6), (17), and (28) 
and the reorder points are translated into days-of-supply using equation (18). 
The processor 303 outputs the computed information to an output device or 
database 304, this information including output days-of-supply for each 
component i€S and output achieved service level for each end product m€M. 

20 Numerical Results 

The purpose of this section is to illustrate numerically the effectiveness 
of the model and the algorithm developed above. 

Problem Data 

Consider a family of desktop computers that consists of six different 
25 end products, denoted PC 1 - PC 6, which are configured from a set of 1 8 



YOR9-2000-0041 



» i 



27 

different components, or building blocks. The components are obtained from 
external suppliers with different supply leadtimes. Table 1 summarizes the data 
used in the numerical study, including the unit cost of each component and the 
supply leadtime (including time for shipment). 

5 Each end product is configured from a set of building blocks as defined 

by its bill-of- material (BOM). The BOM structure is displayed in Table 2. The 
final assembly process, for each end product, takes no more than 1 to 2 hours, 
which is order-of-magnitude shorter than the component leadtime. So this fits 
well with the scenario discussed above. 

10 We remark that the data provided in Tables 1 and 2 are based, on real 

data from the PC assembly system of a major computer manufacturer. To 
protect company confidentiality, we have altered the values of the data. The 
modified values, however, are still relatively close to a realistic product 
representation. 

1 5 We examine the performance of our model under several parameter 

settings. The parameters varied include the budget limit on the expected 
inventory cost, demand variability (in terms of the coefficient of variation), and 
the BOM structure. We also investigate the performance of the model when 
applied in a non-stationary demand environment We check the results obtained 

20 from the optimization against Monte-Carlo simulation. Simulation point 

estimates and confidence intervals were constructed from twenty independent 
replications with 1,000 time periods in each replication. 



Table 1: Unit Costs and Supply Leadtimes of the Components 

Product Cost Component Cost Leadtime 
(in$) (in $) (in days) 

25 PCI 1099 Memory 64MB 290 15 

PC 2 1331 Processor 450 MHZ 307 12 
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Table 2: BOM Structure 



Component Usage 



PCI PC 2 PC 3 PC 4 PC 5 PC 6 
Memory 64MB 111111 
Processor 450 MHZ 1 
Processor 500 MHZ 1 
Processor 550 MHZ 1 1 1 

10 Processor 600 MHZ 1 

Shell 111111 
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Shell (Common Parts 1) 
Shell (Common Parts 2) 
Hard Drive 6.8 GB 
Hard Drive 13.5 GB 
Hard Drive (Common Parts) 
CD ROM 

CD ROM (Common Parts) 
Video Graphics Card 
Ethernet Card 
Software Pre-Load 1 
Software Pre-Load 2 
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1 
1 
1 



1 
1 
1 



1 

1 

1 
1 



1 
1 

1 
1 



1 
1 
1 



Demand Variability 



ki. if 



15 



20 



We first study the effect of demand variability on the service level, 
inventory budget, and optimal base-stock levels. Assume that demand for the 
six end products is independently identically distributed among the products 
and over the time, with a mean of E(D m ) = 100 per period; the coefficient of 
variation is 0.2, 0.4, and 1.0. The budget limits are varied so as to obtain 
service levels (no-stockout probabilities) at about 90%, 95%, and 98%. The 
basic time period in all cases is one day. 

Tables 3 to 5 summarize the numerical results for these cases. The 
column labeled "opt." shows the optimal no-stockout probabilities of the end 
products for t e given inventory budget. The column labeled "sim." shows the 
corresponding 95% confidence interval obtained from simulation. Input to the 
simulation are the base-stock levels for the building blocks generated by the 
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optimization model. The base-stock level of each end product is zero, since 
each end product follows aft configure-to-order operation. 

Table 3: Comparison for End-Product Service Levels; 



E(DJ 100; cv(DJ = 0.2 







B = 800,000 


B = 900,000 


B= 1,000,000 






sim. 


opt. 


sim. 


opt. 


sim. 


opt. 


5 


PCI 


0.949±0.010 


0.923 


0.975±0.006 


0.963 


0.989±0.004 


0.984 




PC 2 


0.949±0.007 


0.923 


0.975±0.005 


0.963 


0.989±0.004 


0.984 




PC 3 


0.980±0.004 


0.923 


0.990±0.003 


0.963 


0.995±0.002 


0.984 




PC 4 


0.966±0.005 


0.923 


0.982±0.004 


0.963 


0.991±0.003 


0.984 




PC 5 


0.991±0.003 


0.945 


0.995±0.003 


0.974 


0.997±0.002 


0.988 


10 


PC 6 


0.945±0.007 


0.923 


0.973±0.005 


0.963 


0.988±0.004 


0.983 



Table 4: Comparison for End-Product Service Levels; 
E(D m )100;cv(2) m ) = 0.4 



B = 1,500,000 B = 1,750,000 B = 2,000,000 





sim. 


opt. 


sim. 


opt. 


sim. 


opt. 


PCI 


0.927±0.013 


0.892 


0.969±0.007 


0.956 


0.989±0.004 


0.984 


PC 2 


0.928±0.008 


0.892 


0.970±0.006 


0.956 


0.989±0.004 


0.984 


PC 3 


0.969±0.006 


0.892 


0.987±0.003 


0.956 


0.995±0.003 


0.984 


PC 4 


0.949±0.007 


0.892 


0.978±0.004 


0.956 


0.990±0.003 


0.984 


PC 5 


0.983±0.004 


0.923 


* 0.993±0.003 


0.969 


0.997±0.002 


0.988 


PC 6 


0.919±0.009 


0.892 


0.967±0.005 


0.955 


0.988±0.004 


0.983 
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Table 5: Comparison for End-Product Service Levels; 
E(DJ 100; cy{DJ = 1.0 



B = 800,000 B = 900,000 B = 1,000,000 





sim. 


opt. 


sim. 


opt. 


sim. 


opt. 


PC 1 


0.926±0.013 


0.923 


0.951±0.011 


0.957 


0.967±0.009 


0.977 


PC 2 


0.926±0.013 


0.923 


0.949±0.008 


0.957 


0.967±0.007 


0.977 


PC 3 


0.954±0.007 


0.923 


0.970±0.005 


0.957 


0.981±0.004 


0.977 


PC 4 


0.933±0.008 


0.923 


0.953±0.007 


0.957 


0.967±0.006 


0.977 


PC 5 


0.971±0.006 


0.945 


0.981±0.005 


0.970 


0.988±0.004 


0.984 


PC 6 


0.928±0.011 


0.923 


0.953±0.009 


0.957 


0.969±0.008 


0.977 



Observe that the no-stockout probabilities returned from the 
optimization model is generally conservative; i.e., lower than the simulated 
values. This is expected as the off-shelf availability of each end product used in 
the optimization model is a lower bound of the true value, as explained above. 
Also observe that the bound becomes tighter as the service level increases, 
along with the increased budget. Furthermore, the bound appears to be closer to 
the simulation results as the coefficient of variation (of demand) increases. 

Among the six end-products, PC 5 is the only one that does not have a 
unique component. Observe that it alone achieves a higher service level than 
the other five end-products, which all have the same service level. In other 
words, all end-products with unique components have binding constraints in 
inequality (24), whereas PC 5 exceeds this service level. 

For the same cases, Table 6 summarizes the comparison for the 
expected overall inventory costs. It is interesting to note that while the 
simulation runs result in higher service levels, the associated inventory costs 
are also slightly higher. (The relative error, however, is still well below 5%.) 
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This suggests that in applications, when the inventory budget is a "hard" 
constraint, we need to input into the optimization model an inventory budget 
limit that is about 5% lower than its real value. 

Table 6: Comparison for Expected Inventory Costs; 
5 E(DJ = 100 





cv(DJ = 0.2 




cv(DJ = 0.4 




cv(DJ = 1.0 




opt. 


sim. 


rel.err. 


opt. 


sim. 


rel.err. 


opt. 


sim. 


rel.err. 


800,000 


841,819 


4.9% 


1,500,000 


1,554,728 


2.5% 


4,000,000 


4,056,800 


1.4% 


900,000 


923,893 


3.5% 


1,750,000 


1,774,716 


1.4% 


4,400,000 


4,419,355 


0.5% 


1,000,000 


1,014,657 


1.4% 


2,000,000 


2,010,165 


0.4% 


4,800,000 


4,792,957 


0.1% 



Table 7 displays the solution found by the optimization model for the 
case of cv(Z) m ) = 1.0. For every component, the optimal solution is shown in 

three forms: base-stock level R* , safety factor k* , and days-of-supply DOS,. 

Clearly, these (optimal) values all increase as, the overall budget increases. The 
1 5 results also illustrate the effect of risk pooling: the higher the degree of 

commonality of a component, the lower the amount of safety stock required to 
maintain the service level. For instance, the safety stock for the 64MB memory 
module, which is used in all six products, varies between 3.86 and 4.48 days- 
of-supply, whereas the safety stock for the ethernet card, which is unique to PC 
20 3, is roughly twice that amount. 
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Table 7: Optimal Solution in Terms of Base-stock Levels, 
Safety Factors, and Days-of-Supply; 
E(2> J = 100; cv(DJ = 1.0 







5 = 4,000,000 




B = 4,400,000 




5 = 4,800,000 






*/ 


DOS, 


R i 


i * 


DOS, 


R i 


K 


DOS, 




Memory 64MB 


11318 


2.44 


3.86 


11509 


2.64 


4.18 


11689 


2.83 


4.48 


5 


Processor 450 MHz 


1760 


1.62 


5.60 


1856 


1.89 


6.55 


1950 


2.16 


7.49 




Processor 500 MHz 


1760 


1.61 


5.59 


1856 


1.89 


6.55 


1950 


2.16 


7.49 




Processor 550 MHz 


4793 


1.99 


3.98 


4940 


2.23 


4.46 


5088 


2.48 


4.96 




Processor 600 MHz 


1745 


1.57 


5.45 


1842 


1.85 


6.42 


1936 


2.12 


7.35 


O 


Shell 


4803 


3.29 


3.00 


4880 


3.43 


3.13 


4940 


3.54 


3.23 


m io 


Shell (Common Parts 1) 


6818 


2.91 


3.36 


6922 


3.06 


3.54 


7045 


3.24 


3.74 




Shell (Common Parts 2) 


6818 


2.91 


3.36 


6922 


3.06 


*3.54 


7045 


3.24 


3.74 


C 


ndiu Lsi ivc u . ovjij 


7251 


2.52 


6 17 


7390 


2.71 


6.63 


7540 


2.91 


7.13 




Hard Drive 13.5GB 


7028 


2.21 


5.42 


7192 


2.44 


5.97 


7356 


2.66 


6.52 




Hard Drive (Common Parts) 


6818 


2.91 


3.36 


6942 


3.09 


3.57 


7045 


3.24 


3.74 


r S 15 


CD ROM 


1863 


2.73 


8.62 


1921 


2.91 


9.21 


1978 


3.09 


9.77 




CD ROM (Common Parts) 


1742 


2.35 


7.42 


1809 


2.56 


8.08 


1877 


2.77 


8.76 




Video Graphics Card 


4538 


2.81 


3.08 


4638 


2.99 


3.27 


4729 


3.16 


3.46 




Ethernet Card 


1742 


2.35 


7.42 


1809 


2.56 


8.0 


1877 


2.77 


8.76 




Software Pre-Load 1 


2740 


2.85 


2.85 


2805 


3.01 


3.01 


2878 


3.19 


3.19 


20 


Software Pre-Load 2 


795 


1.97 


3.94 


842 


2.21 


4.41 


889 


2.44 


4.88 



Unique Components 

Next, we illustrate the special case that each end-product involves at 
least one unique component, in particular, the verification of the optimality 
equations (26) to (29), as discussed at the end of the last section. Recall, in the 
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above example, PC 5 is the only end-product that does not have a unique 
component. To obtain a scenario where each end product has at least one 
unique component, we modify the original bills-of-materials structure by 
letting PC 3 and PC 4 both use the 6.8GB hard drive, while reserving the 
13.5GB hard drive as a unique component for PC 5. This way, PC 5, like all 
the other end-products, also has unique component. 

We run the greedy algorithm on this modified BOM structure fo the 
case of cv(D m ) = 1.0. The results are summarized in Tables 8 and 9. Notice that 
the six end products now achieve the same service level because they all have 
unique components, and thus all constraints in inequality (24) become binding. 

Table 8: Comparison for End-Product Service Levels 
(With Each End-Product Having at Least One Unique Component); 







E(DJ 


= 100; cv(DJ = 


1.0 








B = 4,000,000 




B = 4,400,000 




5 = 4,800,000 






sim. 


opt. 


sim. 


opt. 


sim. 


opt. 


PC 1 


0.934 ±0.012 


0.933 


0.956 ±0.010 


0.963 


0.971 ±0.008 


0.981 


PC 2 


0.932 ± 0.009 


0.933 


0.956 ± 0.007 


0.963 


0.971 ± 0.006 


0.981 


PC 3 


0.961 ± 0.006 


0.933 


0.975 ± 0.005 


0.963 


0.985 ± 0.003 


0.981 


PC 4 


0.940 ± 0.007 


0.933 


0.958 ± 0.006 


0.963 


0.972 ± 0.005 


0.981 


PC 5 


0.946 ± 0.009 


0.933 


0.965 ± 0.008 


0.963 


0.978± 0.007 


0.981 


PC 6 


0.936 ±0.010 


0.933 


0.959 ± 0.009 


0.963 


0.973 ± 0.007 


0.981 



When comparing the results for the original and modified BOM 
structures as shown in Tables 7 and 9, we find that the optimization produces 
similar days-of-supply targets for all components except the two hard drives. 
The safety stock of the 13.5GB hard drive increases from 5.42 to 8.19 days-of- 
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supply, whereas that of the 6.8GB hard drive decreases from 6.17 to 5.05 days- 
of-supply. This again is a consequence of risk-pooling because modifying the 
BOM structure increases the commonality of the high-capacity hard drive, and 
at the same time decreases the commonality of the low-capacity hard drive. 
5 To verify the optimality equations, recall from the discussion above that 

we first determine the Lagrangian multipliers, X m and T|, based on the solutions 
returned by the algorithm, and then check the equations in (26) for all the non- 
unique components. Table 10 lists the gap, labeled "abs.err", between the left 
hand side and the right hand side of the equations in (26) for all the non-unique 
10 components. The gap is indeed very small, on the order of 10" 8 or below. This 
clearly indicates that in this case the greedy algorithm has returned an optimal 
solution. (That there is a gap at all is mainly due to discretization.) The other 
two cases, with the coefficients of variation being 0.2 and 0.4, yield similar 
results. 

1 5 Table 9: Optimal Solution in Terms of Base-Stock Levels, 

Safety Factors, and Days-of-Supply (With Each End-Product 
Having at Least One Unique Component); 
E(DJ = 100; cv(DJ = 1.0 





5 = 4,000,000 




5 = 4,400,000 




B = 4,800,000 






*/ 


K 


DOS 


*; 


K 


DOS 


*/ 


K 


DOS 


Memory 64MB 


11357 


2.48 


3.93 


11542 


2.68 


4.21 


11731 


2.88 


4.55 


Processor 450 MHz 


1781 


1.68 


5.81 


1876 


1.95 


6.75 


1970 


2.22 


7.69 


Processor 500 MHz 


1781 


1.68 


5.81 


1876 


1.95 


6.75 


1970 


2.22 


7.69 


Processor 550 MHz 


4788 


1.98 


3.96 


4936 


2.23 


4.45 


5088 


2.48 


4.96 


Processor 600 MHz 


1767 


1.63 


5.66 


1863 


1.91 


6.62 


1956 


2.18 


7.56 


Shell 


4837 


3.35 


3.06 


4940 


3.54 


3.23 


5038 


3.72 


3.40 
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Shell (Common Parts 1) 


6843 


2.95 


3.40 


6963 


3.12 


3.60 


7080 


3.29 


3.80 


Shell (Common Parts 2) 


6843 


2.95 


3.40 


6963 


3.12 


3.60 


7080 


3.29 


3.80 


Hard Drive 6.8GB 


11525 


2.66 


5.05 


11703 


2.85 


5.40 


11881 


3.04 


5.76 


Hard Drive 13.5GB 


2619 


1.93 


8.19 


2724 


2.18 


9.23 


2827 


2.42 


10.26 


Hard Drive (Common Parts) 


6843 


2.95 


3.40 


6963 


3.12 


3.60 


7080 


3.29 


3.80 


CD ROM 


1860 


2.72 


8.59 


1921 


2.91 


9.21 


1978 


3.09 


9.77 


CD ROM (Common Parts) 


1739 


2.33 


7.38 


1807 


2.55 


8.06 


1877 


2.77 


8.76 


Video Graphics Card 


4561 


2.85 


3.12 


4663 


3.04 


3.33 


4750 


3.19 


3.50 


Ethernet Card 


1739 


2.33 


7.38 


1807 


2.55 


8.06 


1877 


2.77 


8.76 


Software Pre-Load 1 


2752 


2.88 


2.88 


2825 


3.06 


3.06 


2896 


3.24 


3.24 


Software Pre-Load 2 


793 


1.96 


3.92 


842 


2.21 


4.41 


890 


2.45 


4.89 



Table 10: Optimization Results for Non-Unique Components 
(Scenario Where Each End Product Has at Least One Unique 
Component); E(DJ = 100; cv(DJ = 1.0 



B = 4,000,000 B = 4,400,000 B =4,800,000 







* 


abs.err. 


* 

x i 


abs.err. 


* 


abs.err. 


15 


Memory 64MB 


0.0073 


5.20" 8 


0.0041 


3.04" 8 


0.4983 


1.68" 8 




Processor 550 MHZ 


0.0233 


3.02" 8 


0.0129 


1 .84" 8 


0.4949 


1.10' 8 




Shell 


0.0005 


5.59" 8 


0.0003 


3.78' 8 


0.4999 


2.26" 8 




Shell (Common Parts 1) 


0.0018 


5.31- 8 


0.0011 


3.20- 8 


0.4996 


1.73" 8 




Shell (Common Parts 2) 


0.0018 


5.31" 8 


0.0010 


3.20" 8 


0.4996 


1.73" 8 


20 


Hard Drive 6.8GB 


0.0059 


5.21- 8 


0.0034 


3.10- 8 


0.4987 


1.70- 8 




Hard Drive (Common Parts) 


0.0018 


5.31- 8 


0.0011 


3.20" 8 


0.4996 


1.73" 8 




Video Graphics Card 


0.0025 


4.08" 8 


0.0014 


2.58- 8 


0.4994 


1.24" 8 




Software Pre-Load 1 


0.0022 


2.72" 8 


0.0013 


1.76" 8 


0.4995 


0.93' 8 
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Non-Stationary Demand 

We next examine the effectiveness of the algorithm in handling non- 
stationary demand. For instance, when the demand forecast changes weekly, we 
can run the optimization model at the beginning of each week using the 

5 demand offset scheme described above. The base-stock levels are then updated 
weekly, based on the outcome of the optimization. 

To illustrate this approach, consider the product structure described in 
Table L Consider a planning horizon of 100 weeks, over which the mean 
demand for each end-product changes weekly, while the coefficient of variation 

1 0 (which measures the forecast error) remains constant. Specifically, we 

randomly generate the weekly demand from a uniform sample over the interval 
[525,875], while maintaining the coefficient of variation at 2.646. Using the 
algorithm that solves the "dual" problem described in (3 1) and (32), we 
minimize the total expected inventory cost subject to a set of target no-stockout 

1 5 probabilities, 90%, 95%, and 98% (corresponding to p = 0. 1 0, 0.05, and 0.02). 

Table 1 1 compares the achieved service levels and the total expected inventory 
cost between simulation and the optimization model. Input to the simulation 
are the component base-stock levels, by week, generated by the optimization 
algorithm. 

20 Table 11: Comparisons for End-product Service Levels and 



Total Inventory Cost in Non-Stationary Demand Environment 





P = 


0.10 


P = 


0.05 


P = 


0.02 


sim. 


opt. 


sim. 


opt. 


sim. 


opt. 


PC 1 


0.928±0.009 


0.900 


0.952±0.007 


0.950 


0.971±0.004 


0.980 


PC 2 


0.926±0.007 


0.900 


0.949±0.005 


0.950 


0.969±0.003 


0.980 


PC 3 


0.944±0.006 


0.900 


0.962±0.004 


0.950 


0.975±0.003 


0.980 
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PC 4 


0.929±0.007 


0.900 


0.951±0.0Q6 


0.950 


A AiCO_i_A (\C\A 

O.VOo±U.UU4 




PC 5 


0.950±0.006 


0.900 


0.965±0.004 


0.950 


0.978±0.003 


0.980 


PC 6 


0.922±0.009 


0.900 


0.950±0.008 


0.950 


0.970±0.005 


0.980 


Total cost 


$4,039,913 


$3,830,505 


$4,505,698 


$4,329,152 


$5,077,240 


$4,924,748 



5 For 90% and 95% service level targets, the simulated no-stockout 

probabilities do exceed the values returned by the optimization. At 98% 
service level, the simulated values are slightly below their analytical 
counterparts. This is largely due to the weekly update of the base-stock levels: 
whenever there is an increase in demand, and hence an increase in base-stock 
l ;f 10 levels, from one week to the next, there is a time lag in building up the 

VI inventory, due to the supply leadtime, and the available stock may be 

12 insufficient to accommodate the increase in demand. This shortfall can be 

eliminated if the base-stock levels are updated more frequently, say, daily. We 
1 0 further observe that the analytical estimates of the expected total inventory 

1% 15 cost are within 5.5% of the simulated values. 

* z Figure 4 shows the optimal safety-stock levels (i.e., base-stock minus 

M work-in-progress (WIP)) converted to days-of-supply (DOS), for a subset of 

components in the example. Notice that the DOS values fluctuate very little 
over time. This suggests, in situations where the forecast error (in terms of its 
20 coefficient of variation) remains relatively constant over time, it suffices to 
follow a simple policy that tracks, for each component inventory, a constant 
DOS target (instead of updating the base-stock levels, which keep changing 
from period to period). Further, the results clearly demonstrate the benefit of 
risk-pooling: building blocks that go into a larger number of end products 
25 (e.g., memory modules) require much less safety stock than building blocks 

that are only used in a small number of end products (e.g., CD ROMs). 
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Process Re-engineering Application 



Here we describe the study mentioned earlier, which was part of a 
larger project aimed at the re-engineering of an existing business process from 
a build-to-stock operation to an assembly-to-order operation centered around 
"building blocks" (i.e., keeping inventory only at the component level). 
To carry out the study, we have developed two basic models: the "as-is" 
model that is a reflection of the present operation, and the "to-be" model that 
is based on the optimization model described in the previous sections - in 
particular, with the component inventory levels generated by the algorithm. 
For both models, we aggregate the business' production-inventory system into 
two stages, the first stage consists of the components, or building blocks, and 
the second stage includes the assembly and the order fulfillment. 

Three factors have been identified as focal points of our study: 

(i) manufacturing strategy - the "as-is" operation versus the "to-be" 
model; 

(ii) the accuracy of demand forecast at the end-product level versus at the 
component level; and 

(iii) the effect of mass customization as a result of, for instance, direct sales 
over the internet 

To study the first factor, we select a major product family at the 
business, which consists of 18 end-products configured from a total of 17 
components. We use the business' existing data - including BOM, unit costs, 
assembly and procurement leadtimes - to run a detailed simulation model. The 
demand for each end-product is statistically generated, based on historical 
data. The days-of-supply targets are set to meet a required service level of 
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95% (for all end-products). Following the business' current practice, these 
targets are determined using a combination of simple heuristic rules and 
judgement calls from product managers, and verified by simulation (via trial 
and error). 

5 We then feed the same data, including the same statitically generated 

demand streams, into the optimization model, which eliminates the finished- 
goods inventory at the end-product level and optimally sets the base-stock 
level for each component inventory. The optimization model minimizes the 
overall inventory cost while meeting the same service level of 95% for all end- 

1 0 products. We take the optimal base-stock levels, and rerun the simulation to 
verify the analytical results. 

Figure 5 shows the comparison between the "as-is" and the "to-be" 
models, in terms of the overall inventory cost. (To protect proprietary 
information, the vertical axis in all figures is normalized with respect to the 

15 inventory cost of the "as-is" model, which is 100.) As expected, the inventory 
cost at the end-product level is virtually eliminated in the "to-be" model. (The 
cost shown is due to WIP; the cost due to finished goods is nil) In contrast, 
the "as-is" model keeps a significant amount of end-product inventory. On the 
other hand, the amount of component inventory is higher in the "to-be" model, 

20 which is again expected, since the required service level of 95% is common to 
both models. Overall, the "to-be" model reduces the overall inventory cost by 
about 30%. 

Note in the above study, both models use the same demand forecast, at 
the end-product level. The "to-be" model, however, can easily switch to 
25 forecasting demand directly at the component level. This will result in 

improved forecast accuracy, as we can take advantage of parts commonality, 
as each component is generally used in several end products. Hence, in our 
study of the second factor, we evaluate the effect of forecast accuracy through 
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a sensitivity analysis. Figure 6 shows the overall inventory cost associated 
with three different levels of forecast accuracy. The first two columns, repeat 
the comparison in the last figure; i.e., both ("as-is" and "to-be" models assume 
30% forecast error (i.e., the coefficient of variation equals 0.3) at the end- 
product level; the next two columns represent improved forecast errors, at 
20% and 30%, achieved by the "to-be" model through forecasting at the 
component level. 

Our study of the third factor aims at analyzing the impact on inventory 
when the system supports a richer product set, in terms of product variety. The 
motivation is to support mass customization. In the Internet-based, direct-sales 
environment, for instance, the number of different product configurations that 
customers want to order can be significantly larger than what is currently 
supported in the build-to-stock environment. Figure 7 shows the inventory 
costs: the four columns on the left correspond to the current product set (lx), 
with the first scenario (SI) being the "as-is" model, and the other three being 
the "to-be" model at the current (S2) and improved (S3, S4) forecast accuracy 
levels, respectively; the four columns on the right repeat these scenarios with a 
product set that is ten times as large in variety (lOx), while maintaining the 
overall volume. (Also refer to Table 12 for a summary of all the different 
scenarios.) 

Observe that as the product variety increases, a significantly higher 
level of inventory is required in the "as-is" model. This is because forecast 
accuracy will deteriorate when the end products proliferate (i.e., larger 
varieties at smaller volumes). On the other hand, in the "to-be" environment, 
the increase in inventory cost is very modest. This is because the proliferation 
of end-products will have minimal effect on the forecast accuracy at the 
building-block level, due to parts commonality. This strongly supports the fact 
that the building-block model is the right process to support an internet-based, 
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direct-sales operation. 



Table 12: Summary of the Scenarios Used to 
Study the Effect of Product Variety 



Scenario 


Description 


lx Cases 


lOx Cases 


SI 


"as-is" 


original product set, 
jU/o iorecdbL crrur, 
90% service 


ten times larger product set, 
30%xy f 10 forecast error at 
MTM level 


S2 


"to-be" 


forecast at MTM 
level, 30% forecast 
error, 90% service 


ten times larger product set, 
forecast error as in S2(lx) 


S3 


"to-be" 


forecast at BB level, 
20% forecast error, 
90% service 


ten times larger product set, 
forecast error as in S3(lx) 


S4 


"to-be" 


forecast at BB level, 
10% forecast error, 
90% service 


ten times larger product set, 
forecast error as in S4(lx) 



The examples given demonstrate the advantages of the invention 
1 0 applied to the business of manufacturing personal computers (PCs). However, 
the invention has applications to other configure-to-order (CTO) 
manufacturing systems and is not limited to manufacturing PCs. Thus, while 
the invention has been described in terms of a single preferred embodiment, 
those skilled in the art will recognize that the invention can be practiced with 
1 5 modification within the spirit and scope of the appended claims. 
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CLAIMS 



Having thus described our invention, what we claim as new and desire 
to secure by Letters Patent is as follows: 



1 1 . A method of managing manufacturing logistics of end products comprising 

2 the steps of: 

3 maintaining an inventory of components, which components, termed 

4 "building blocks", are built to stock; 

5 configuring-to-order end products using said components; and 

6 replenishing said components from suppliers following a base-stock 

7 policy. 

1 2. The method of managing manufacturing logistics of end products recited in 

2 claim 1 , wherein the end products are personal computers (PCs) and the 

3 components are stock computer components. 

1 3. The method of managing manufacturing logistics of end products recited in 

2 claim 1 , wherein the base-stock levels are derived from a greedy algorithm 

3 which iteratively reduces inventory budget until a budget constraint is 

4 satisfied. 

1 4. A computer implemented process of managing manufacturing logistics of 

2 configure-to-order end products comprising the steps of: 

3 a) initializing the process by setting x t := 0 for each / E S, setting i mi := 

4 ?(X mt > 0), setting P m := 0 for each m6M, and setting P := 0, where S is a set 

5 of components indexed by /, Mis a set of end products indexed by m, x, is the 

6 probability of no-stockout of a component of index /, r m is the probability 
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7 that a positive number of units of component / is used in the assembly of an 

8 end product indexed by m, p m is the probability of stockout of an end product 

9 of index m, and P is the upper limit on the stockout probability over all end 

10 products; 

11 b) setting the set of active components to A :== { } ; 

12 c) considering each i E S, followed by considering each end product m 

13 that uses component i in its bill-of-material; 

14 d) setting p w := p m + r m A, for all m such that / E S m where A is a unit 

15 step size; 

16 e) computing the difference S, := max /7I {p ;w } - p; 

17 f) determining if S, <> 0, and if so, then adding component index i to the 

1 8 set of active components, A := A + {*} ; 

19 g) determining if the set of active components is non-empty, and if so, 

20 then setting B := A, otherwise setting B := S where B is a set of component 

21 indexes; 

22 h) finding z* := arg max, ei? {-c,a y g'(x, + A/2)}, where -#'(•) follows the 

- t _ i \- x 

23 equation -g'Qe) ~ -$(0 (x))- — = = — = , where $(.) is the 

24 probability distribution function of the standard normal variate, and is 

25 the probability density function of the standard normal variate; 

26 i) setting x,* := + A to update the probability of no-stockout of 

27 component/*; 

28 j) computing P := max mej J3 m , and checking whether inequality 

29 ]T) c Pfg( x ) * By where B is the budget limit on the expected overall 

30 inventory cost, is satisfied and if so, stop; 

31 k) otherwise, updating p m and for each m E M t *, set p m := P m + r mi A, and 
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going to step b). 

5. A system for managing manufacturing logistics of end products comprising: 

means for maintaining an inventory of components, which 
components, termed "building blocks", are built to stock; 

means for configuring-to-order end products using said components; 

and 

means for replenishing said components from suppliers following a 
base-stock policy. 

6. The system for managing manufacturing logistics of end products recited in 
claim 5, wherein the end products are personal computers (PCs) and the 
components are stock computer components. 

7. The system for managing manufacturing logistics of end products recited in 
claim 5, wherein the base-stock levels are derived from a greedy algorithm 
which is iteratively computed by a processing unit to reduce inventory budget 
until a budget constraint is satisfied. 

8. A method that translates end-product demand forecast in an 
assemble-to-order (ATO) environment into a forecast for components, taking 
into account outbound leadtime comprising the steps of: 

defining the demand D m (t) of type m in period /, each unit of type m 
demand requiring a subset of components, denoted S m £ S, as 

m€M l 
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7 deriving mean and variance for demand £>,(/) as 

8 V[Dfffl = E EE[^+0]P[C-«L 311(1 
Var[D£)] = E EE[Z)if + {)]P[C^] 

-E(EEP»(M)]P[xr^]) 2>Kspeoovely - 

1 9. The method recited in claim 8, wherein the ATO environment is extended 

2 to a configure-to-order (CTO) environment for stationary demand, taking into 

3 account batch sizes comprising the steps of: 

4 translating end-product demand into demand for each component / (per 

5 period) as 

7 where X m (k) 9 for k = 1, 2,..., are independent, identically distributed (i.i.d.) 

8 copies ofX mi ; 

9 deriving marginal distributions, and then the mean and the variance of 

10 X m ZB 

11 E[Z>J = E E[XJE[Z)J, and 
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VarPJ = E (EPJVar[^.] + Var[i>JE 2 [Zj) 

12 m6jW/ / 

= E E 2 [Z JE[2)2l + VarEXJE[Z) J-E 2 [X ,.]E 2 [Z> J> res ? ectlvey - 



1 10. The method recited in claim 9, extended to non-stationary demand, 

2 wherein the mean and the variance of X ml are generalized as 

3 W>m = E E[XJEeP ot (^0]P[C=«] 5 and 



o 4 

[f! Var[D,(0] = E E 2 (XjE E[Z>>{)]P[C t= «] 

P 5 + E Var(^ f )E E[D w (r + «)]P[C 1 ^] 

!: - E E 2 tUf E E[D m (^{)]P[Z,: ut =C]) 2 , respectively. 

s* "Hi 

r — 



1 11. The method recited in claim 9 ? further comprising the steps of: 

2 defining R t (t) as a reorder point (or, base-stock level) in period / as 

3 R t (t) := ji f (0 + «t)O f (0» 

4 where k t (f) is the desired safety factor, while and G t (f) can be derived (via 

5 queuing analysis as 
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li.(t) = £ E[Z>^)], and 



a f(0 = z2 Varfi).^)], respectively, 



where 4 I . := E[L i ] is expected in-bound leadtime; and 

translating R t (t) into "days of supply" (DOS), where the ^(t) part of 
i? r (0 translates into periods of demand and the kfflofa) part of R t (t) is turned 
into 



periods of demand so that R t (f) is expressed in terms of periods of DOS as 



dos/o = r, 



in 



off) 



12. The method recited in claim 11, wherein demand is stationary in which for 
each demand class m, DJJ) is invariant in distribution over time, so that the 
mean and the variance of demand per period for each component i reduce to 
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= {™E[£>.] 5 and of = { "VarfDJ , respectively, and 



R ( = l*E[Dyk f ffi sd[iy , and hence, 



R, 



DOS,. = — i- = C+kfl t Jl? = F 
1 E[D t ] ' 1 ' V ' 



7 
8 



where 0, := sd[DJ/E[DJ is the coefficient of variation of the demand per 



period for component i, and hence 8 . A/{ j™ is the coefficient of variation of 



the demand over the leadtimeC" 1 . 



1 
2 
3 
4 
5 
6 
7 
8 



13. A method that relates service requirements to base-stock levels of the 
components in an assemble-to-order (ATO) environment comprising the steps 
of: 

defining each order of type m as requiring exactly one unit of 
component i € S m , a as a required service level, referred to as off-shelf 
availability of all the components required to configure a unit of type m 
product, for any m, and E, as an event that component / is out of stock; 

determining a probability P for each end product mE M, 



P[u,- S £ ( .] s l-o, and 



10 



n" ieS E,l = T,^E)-Y t ^iEpE)^V(EpEpE k )- • • •, and 
* i i<j i<j<k 
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14. The method recited in 13, wherein the method is extended to a configure- 
to-order (CTO) environment taking into account batch sizes, further 
comprising the steps of: 

defining A c S m as a certain configuration, which occurs in a demand 
stream with probability ?(A); 

weighting a no-stockout probability, Y[ i€A $0O> by P(4); 

changing the service requirement to 



a s £ ?(A) n*(*!) 



- EpW [l-E *(*»)] 



AcS m ieA 




extending the CTO environment the service requirement to 



15. A method that translates service requirements in terms of leadtimes 
into requirements for the off-shelf availability of the components comprising 
the steps of: 
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relating an off-shelf availability requirement to standard customer 
service requirements expressed in terms of leadtimes, W m , where a required 
service level of type m demand is 



nw m zwj z a, meM, 



where w m 's are given data and P is probability; 

when there is no stockout at any store i 6 S m , denoting the associated 

probability as X Qm (f), a delay being Z?*, the out-bound leadtime; 

when there is a stockout at one or several stores in the subset s c S m9 
denoting the associated probability as Ks m (f), so that the delay becomes 

£° Ut + t*> where T s is the additional delay before the missing components in s 
become available; 

determining V[W m swJ = ^PtC^J + E * JW^+Vi* 



m 

nx 



and 

assuming that 



r OUt j T Out — ^ 

L m <, w and L m +t >w 



both hold almost surely, so that when the (nominal) outbound leadtime is 
nearly deterministic and shorter than what customers require, whereas the 
replenish leadtime for any component is substantially longer. 
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LARGE INVENTORY-SERVICE OPTIMIZATION IN 
CONFIGURE-TO-ORDER SYSTEMS 



ABSTRACT OF THE DISCLOSURE 



A manufacturing process is migrated from an existing operation to a 
configure-to-order (CTO) system. As the CTO operation will eliminate the 
"machine-type model" (MTM) inventory of the existing operation, the 
emphasis is shifted to the components, or "building blocks", which will still 
follow the build-to-stock scheme, due to their long leadtimes, and hence still 
require inventory. The solution involves an inventory-service trade-off of the 
new CTO system, resulting in performance gains, in terms of reduced 
inventory cost and increased service level. Other benefits include better 
forecast accuracy through parts commonality and risk-pooling, and increased 
customer demand, as orders will no longer be confined within a restricted set 
of pre-configured MTMs. 



Components End products 

S={1,2,...,7} M={1,2,3} 




S, ={1,2,3} 
S 2 = {2,4,5} 

S 2 = {3,5,6,7} 



Figure 1 
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Figure 3B 
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Docket No.; YOR9-2000-0041 

Application for United States Patent 
Declaration and Power of Attorney 

As a below named inventor, I hereby declare thai: 

My residence, post office address and citizenship are as stated below next to my name; 

I believe I am an original, first and joint inventor of the subject matter which is claimed and for which a patent 
is sought on the invention entitle d LARGE TNVKNTQ R Y-SERVICE.OPTTMTZATION IN CONFIGURE-TO- 
ORDER SYSTEMS t he specification of which: 

(check & is attached hereto 

one) 

n was filed on as 

Application Serial No- 

and was amended en (if applicable) 

I hereby state that I have reviewed and understand the contents of the above identified specification, including the 
claim?, as amended by any amendment referred to above. 

-~ I acknowledge the duty to disclose information which is material to the examination of this application in accordance 

Hfith Title 37, Code of Federal Regulations, § 1.56(a) * 

W I hereby claim foreign priority benefits under Title 35, United States Code, §119 of any foreign applications) for patent 
IMr inventor s certificate listed below and have also identified below any foreign application for pateni or inventor's certificate 
^having a filing date before that of the application on which priority is claimed: 

Lftiar Foreign Application^) Priority Claimed 



s CNumber) (Country) (Day/Monrh/Year Filed) yes no 

*(Nomber) (Country) ~~ (Day/Month/Year Filed) yes no 

I hereby claim the benefit under Tide 35, United States Code, $ 120 of any United States appiieation(s) listed below and, 
insofar as the subject matter of each of the claims of this application is not disclosed in the prior United States application in the 
fanner provided by the first paragraph of Title 35, United States Code, § 1 12, 1 acknowledge the duty to disclose material 
-information as defined in Title 37, Code of Federal Regulations, § 1 .56(a) which occurred between the filing date of the prior 
application and the national or PCT international filing date of this application: 



(Application Serial No,) (Filing Date) (Status; patented, pending, abandoned) 

Power of Attorney: As a named inventor, I hereby appoint Manny W. Schecter, Reg, No. 31.722, Terry J. Ilardi. Reg, 
No. 29,936, Stephen C. Kaufman, Reg. No. 29,531, Louis J. Pereello, Reg. No. 33,206, Jay P. Sbroilini, Reg. No. 36,266, Robert 
M- Trcpp, Reg. No. 25,933, Daniel P. Monis. Reg. No. 32.053, Wayne L. Ellenbogen. Reg. No. 43,602. Douglas W. Cameron. 
Reg. No. 31,596, David M, Shofu Reg. No, 39.835, Christopher A- Hughes. Reg. No. 26,914, Bdward A. Pennington, Reg. No. 
32,588, John E, HoeL Reg. No. 26>279- Joseph C. Redmond, Jr., Reg. No 1S.753, C. Lamont whitham. Reg. No. 22,424 7 
Marshall M. Curtis, Reg, No, 33,138, and Michael E. Whitham, Reg, No. 32,635, as attorneys and/or agents to prosecute this 
application and transact all business in the Patent and Trademark Office connected therewith. All correspondence should be 
directed to Whitham, Curtis & Whitham, Reston International Center, 1 1800 Sunrise Valley Drive, Suite 900, Reslon, Virginia 
20191. Phone calls should be directed to Whitham, Curtis & Whitham, at 703/39 J -25 10, 
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I hereby declare dim oil statements made herein of my awn knowledge ere true end that all eutfMOttots msdtr on information end 
belief at* be1i6ved ro be true; and farther thai ihese Biaicmwib wcie nude with th» knowledge tku willful fako sinioments and dm 
like 44 (hide bttf punishable by fift* Of Imprisonment, or boih) Oflder Boouon l6M of Tide 1 B of the Untied SuU* Code and lhp( 
such wiltftil fake sLucmonu may jeopardize die validity of the eppUearion or any patem issued thereon. 

(1) Inventor: 




Signitw* r ^ ^ l^U— V- D«e: 

'RmMmwt; 19 PoAduliw Ctoea. Chappiqiu, NY 1 
Cniwnshlp; Cututa 

Pot( Qfffco AMmt; $WW •* Retidenoe 

(2) Inventor. MarttsEnl G~~YIa I 

Signature; Y\yzkK*~\ km _ white* 

Residence: 1B70 Baldwin Ro*d AS, Yorloown Ifelghu, MY IQ59& 

CitucnsMp: Germany 

Post Office Address: Some u Residence 

(3) Inventor; Qrace Yuh-Jiuo Lin 

Rosfdcneof P Oorey J>rive> Chappeque. NY 105 14 
Citizenship; Taiwan. ROC 
1*011 Office Address: Same as Reeldonec 
W Inventor; Psvtd Ps*We| yso 

Resideecc; 1261 Underbill Avenue, Vorktown Heights, NY 10S9S 
ChJaaaihip: Unload Suae* of America 
Cost OfKoe Address; Sumo es Resideoeo 
■TWo 37, Code of Federal Rc^ukdons, Jl-56<a); 

M A dmy of candor end good feidi toward too Front ami Trademark Office rests oo die inventor, on caeh attorney or teem who 
prepares or prosecutes the iippSiciidon end on every other individual who is substantively Evolved in the preparation or 
pwncwJw, w f ih. w lloatlo W «nd who fa iissoeliitod wl A J* l nyonu?rf w |d, ^ ttV ^nee or wldj anyone lu whom di«e b on 
0 bb^UoA 10 ^ w Uc«lo«. All W h indmdwds Iwe a doty to disclose lo 1U0 Offic* iMfounedo-i ibey are ewero of 
vft.cn « ma«eri 0 | to the esamJfladon of iho «pptieeiioii. 4udi ioforo^on k n^cdat where ifcere 1» substantial likelihood that a 
rccaonible examiner *o«W consider U imports {^0^^,- whether to Uic *pp««lioa to iwuc a s a paieou The*»y is 
commenAOfftie with the degree of mirolYcment in tho proparnuon cr pto^eution of the appHcation. 
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(b) Under this section information is material to patentability wlien it is not cumulative to information already of record or being 
made of record in the application, and (1) it establishes, by itself or in combination with other information, a prima facie case of 
unpatentability; or (2) it refutes, or is inconsistent with, a position the applicant takes in: (i) opposing an argument of 
unpatentability relied on by the Office, or (ii) asserting an argument of patentability. 
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